Perfect and near perfect adaptation in a model of bacterial chemotaxis 



(N 
O 

o 

(N 

> 
O 

o 

(N 



X3 

Of 



Bernardo A. Mello 1 - 2 and Yuhai Tu 1 

'IBM T.J. Watson Research Center, P.O. Box 218, Yorktown Heights, NY 10598 
A Physics Department, Catholic University of Brasilia, 72030-170, Brasilia, DF, Brazil 

(Dated: February 2, 2008) 

ABSTRACT 

The signaling apparatus mediating bacterial chemotaxis can adapt to a wide range of persistent 
external stimuli. In many cases, the bacterial activity returns to its pre-stimulus level exactly and 
this "perfect adaptability" is robust against variations in various chemotaxis protein concentrations. 
We model the bacterial chemotaxis signaling pathway, from ligand binding to CheY phosphorylation. 
By solving the steady-state equations of the model analytically, we derive a full set of conditions 
for the system to achieve perfect adaptation. The conditions related to the phosphorylation part 
of the pathway are discovered for the first time, while other conditions are generalization of the 
ones found in previous works. Sensitivity of the perfect adaptation is evaluated by perturbing 
these conditions. We find that, even in the absence of some of the perfect adaptation conditions, 
adaptation can be achieved with near perfect precision as a result of the separation of scales in 
both chemotaxis protein concentrations and reaction rates, or specific properties of the receptor 
distribution in different methylation states. Since near perfect adaptation can be found in much 
larger regions of the parameter space than that defined by the perfect adaptation conditions, their 
existence is essential to understand robustness in bacterial chemotaxis. 
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I. INTRODUCTION 

The motion of coliform bacteria (such as E-coli) is 
driven by rotation of several fLagella attached to the 
cell body. When the flagella rotate counter-clockwise 
(CCW), the flagella form a bundle that pushes the bac- 
terium in a smooth motion (runs) with a high degree of 
directionality. On the other hand, when the flagella ro- 
tate clockwise (CW), the flagella bundle flies apart and 
the bacterium tumbles, randomizing the direction of the 
subsequent run. The frequency with which the tumbling 
motion occurs decreases with increasing concentration of 
attractant (or decreasing concentration of repellent). As 
the result, the bacterium performs a biased random walk 
towards higher concentration of attractant. This mecha- 
nism gives the bacterium its ability to follow the gradient 
of chemical concentration, i.e., chemotaxis. 

From the sensing of external stimulus to the activation 
of motor regulator protein, a series of chemical reactions 
are involved in relaying and regulating the signal. For 
recent review on bacterial chemotaxis signaling pathway, 
see references (Falke, 1997; Bren, 2000; Bourret, 2002). 
The major players in the chemotaxis signal transduction 
pathway are the transmembrane chemotaxis receptors 
and 6 cytosolic proteins: CheA, CheB, CheR, CheW, 
CheY and CheZ. The receptor forms a complex with 
the histidine kinase CheA through the adaptor protein 
CheW. The receptor has a ligand binding domain located 
at the periplasm to sense the external signal, such as the 
concentration of attractant (or repellent). The activity 
of CheA is affected by the properties of the receptor, for 
example, whether the receptor is ligand bound or not. 
When chemoattractant binds to receptor, CheA activ- 
ity is suppressed. The histidine kinase CheA, once acti- 
vated, acquires a phosphate group through autophospho- 
rylation, and subsequently transfers the phosphate group 



to the response regulator protein CheY or the demethyla- 
tion enzyme CheB. The phosphorylated CheY (CheY-P) 
then interacts with the motor and increases the motor's 
CW rotation bias. This is the "linear" signal transfer- 
part of the bacterial chemotaxis pathway. Like many 
other biological sensory systems, the bacterial chemo- 
taxis pathway also has the ability to adapt to persistent 
external stimulus. The adaptation in bacterial chemo- 
taxis is facilitated by the methylation and demethylation 
of the receptor, which serves as the feedback control of 
the system. The methylation and demethylation pro- 
cesses are catalyzed by CheR and ChcB-P respectively 
and are slow in comparison with the other reactions. 

Because of the excellent understanding of each individ- 
ual reaction of the pathway, mathematical modelling of 
bacterial chemotaxis signal transduction has been very 
fruitful (Bray, 1993; Hauri, 1995; Barkai, 1997; Spiro, 
1997; Morton-Firth, 1998; Morton-Firth, 1999; Yi, 2000). 
Besides being useful in understanding specific aspects of 
chemotaxis experiments, modelling is essential in gain- 
ing insight about general properties of biochemical net- 
works. One important general problem is to understand 
the functional stability of biochemical networks under 
changes of various pathway parameters, such as concen- 
trations of enzymes and reaction rates. Parameter fluc- 
tuations are inherent for biological systems in the real 
world, so robustness, i.e., the insensitivity of important 
system properties with respect to parameter variation 
and fluctuation of protein concentrations, is crucial for 
the proper functioning of the biological systems. 

Experimentally, it was observed that after initial re- 
sponse to some external stimulus, such as sudden changes 
of aspartate concentrations, the bacteria tumbling fre- 
quency often reverts to its original value with high ac- 
curacy, independent of the strength of the external stim- 
ulus (Berg, 1972). This accurate adaptation is gener- 
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ally believed to contribute to the high sensitivity of bac- 
terial chemotaxis to a wide range of external stimulus 
(5 orders of magnitude). In a recent work, Barkai and 
Lciblcr (Barkai, 1997) investigated the robustness of per- 
fect adaptation in bacterial chemotaxis, they used a two- 
state (active or inactive) model (Asukura, 1984) for the 
receptor complex in explaining the phenomena. In their 
model, they assumed that CheB only demethylates active 
receptors, whereas CheR methylates all receptors indis- 
criminately. They showed, by extensive simulation of the 
two-state model, that as long as the above conditions are 
satisfied, adaptation is achieved with high precision, inde- 
pendent of specific values of the rate constants or enzyme 
concentrations. In a subsequent study, Alon et al. (1999) 
provided experimental evidence for the robustness of the 
perfect adaptation over large variations in chemotactic 
protein concentrations. 

The Barkai-Leibler (BL) model clearly captured one of 
the essential ingredients for perfect adaptation in bacte- 
rial chemotaxis. Recently, Yi, Huang, Simon and Doyle 
(Yi, 2000) further studied the Barkai-Leibler model ana- 
lytically, and summarized all the conditions for perfect 
adaptation within the BL model beyond those identi- 
fied in the original paper. However, the BL model is 
a simplified description of the real chemotaxis pathway. 
For example, the BL model neglects the phosphorylation 
part of the pathway altogether and assumes the satura- 
tion of methylation enzyme CheR, which is questionable 
(Morton-Firth, 1999). 

In this work, we study a more complete model of the 
chemotaxis signal transduction pathway, similar to the 
deterministic version of the model proposed by Morton- 
Firth et al. (1998), where both the methylation and 
phosphorylation processes are taken into account. Our 
goals are to understand whether (mathematically) per- 
fect adaptation, defined as when steady-state CheY-P 
concentration is independent of ligand concentration, can 
be achieved for the full model, and to identify the condi- 
tions for such perfect adaptation. The sensitivity of the 
perfect adaptability, or robustness, is then studied by 
perturbing these conditions. Such study can help us un- 
derstand adaptation in real biological systems where not 
all the perfect adaptation conditions are satisfied, it can 
also provide possible explanations for cases where per- 
fect adaptation is not achieved, e. g., for serine response 
(Berg, 1972). 



II. MODEL 

For the purpose of this study, we consider only those 
receptors that form complex with CheW and CheA. We 
label the receptor complex by T n \, where n(G [0,4]) is 
the number of methyl groups added to the receptor and 
A (= o, v) represents the ligand occupied (o) and vacant 
(v) state of the receptor. Superscripts are also used to 
describe whether the receptor complex is phosphorylated 
(P) or un-phosphorylated (U), bound to CheR/CheB-P 



TABLE I: Chemical species and subspecies. Total concentra- 
tions are taken from reference (Morton-Firth, 1999). 



Species Description 



Concentration 



[T 



[R 1 



[B 1 



[TnX] 

[T F ] 
[T P ] 
[T U ] 

[R F ] 

[B F ] 
[B P ] 
[B PF ] 

[Y P ] 



Total taxis aspartate receptor (Tar) 2.5 fiM 

Receptor with n methyl groups, ligand 

binding site occupied (A = o) or vacant (A = v) 

Free (CheR and CheB unbound) receptor 

Phosphorylated receptor 

Unphosphorylated receptor 

CheR 0.176 ^M 

Free (not bound to T) CheR 

CheB 2.27 

Free (not bound to T) CheB 

Phosphorylated CheB 

Free phosphorylated CheB 

CheY 18 

Phosphorylated CheY 



TABLE II: Chemotaxis signal transduction reactions 



Ligand binding T nv + L <-> T n L(= T no ) 



Methylation 



T n + R F <- 
T n + B PF 



-» T n B p 



T„R — 
T„B P 



Tn+l + 

-> T n _i 



R 



F 
B PF 



Phosphorylation -* T£ 
T P + Y U 



T" + Y p Y 
T p + b uf -» + B PF B 1 



Y u 

-> B UF 



or free (F). Superscript (T) is used to label total con- 
centrations of different proteins. The superscripts are 
not mutually exclusive, e. g., [B PF ] is the concentration 
of phosphorylated free (not bound to receptor) CheB. In 
table |, some of the chemical species of the chemotaxis 
pathway are shown, where the values of the total con- 
centrations are taken from (Morton-Firth, 1999), except 
for the total CheR concentration, which we have reduced 
slightly in order to have the same average methylation 
level as reported in (Morton-Firth, 1999), where recep- 
tors other than Tar were included in the simulation. 

The bacterial chemotaxis pathway can be divided into 
3 processes: receptor ligand binding, receptor methyla- 
tion/demethylation and phosphorylation of CheA, CheB 
and CheY. The reactions involved in each of the three 
processes are listed in table [lj|. Since the ligand bind- 
ing process is much faster than the other two, the lig- 
and binding reaction can be considered to be always in 
quasi-equilibrium. The receptor's ligand binding status 
directly affects both the CheA auto-phosphorylation rate 
and the receptor methylation/demethylation rates. The 
CheA auto-phosphorylation rate is also affected by the 
methylation state of the receptor. Finally, since only the 
phosphorylated CheB can efficiently demethylate the re- 
ceptor, the methylation process is also affected by the 



FIG. 1: Illustration of the methylation and demethylation 
reaction network, n is the methylation level of the receptor . 



an open question, it is supported by some exper- 
iments (Shapiro, 1994; Shapiro, 1995). The net- 
work of methylation/dcmcthylation reactions are 
illustrated in Fig. |[ 

If we assume the above reactions follow Michaelis- 
Menten kinetics and the dissociation rates for the 
bound state are independent of A, i.e., whether the 
receptor is ligand bound or not, the bound state 
concentration can be written as: 



where K% = [(1 - L n )K* v 1 + L n Kg Y l is thc 



Michaelis constant of the combined (vacant and 
ligand-bound) receptor state and the superscript 
F denotes the free enzyme and the free substrate 
(receptor) concentrations. 

Since the receptors and the enzymes can exist either 
in their free form or bound to each other, the total 
concentrations of enzymes, and the concentration 
of receptors with n methylated sites are given by 
the following equations: 



[T n E] 



[T r f}[E F ] 



(6) 



[B P ] 




n=0 n 




(7) 



[T n 



1 + 



[R F ] , [B PF )\ 



K R K B 1 



where the [R T ], [B p ] and [T n ] are the concentra- 
tions of CheR, phosphorylated CheB and receptors 
with n methyl groups, respectively. 
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The kinetic equation for the receptor concentra- 
tions [T n ] at each methylation level can be written 
as: 



d[T n ] 
dt 



J n i 



(10) 



where J n is the net flux from methylation level n to 
level (n + 1), which is just the difference of methy- 
lation and demethylation rates between these two 
states. Using the bound state concentration given 
in Eq. [| J n can be written as : 



Jn 



\R F Wi 



[B PF }[T F +1 ] 
K* +1 



(11) 



(0 < n < 3), 



where fc„ and k B are the catalytic constants for 
the methylation and demethylation reaction respec- 
tively, which are assumed to be independent of 
A, the ligand binding status of the receptor. The 
boundary conditions for the methylation flux are: 
J_i = J 4 = 0. 

3. The auto-phosphorylation of CheA reaction is given 
by: 



T 



(12) 



the phosphate group is subsequently transferred 
from CheA-P to CheB and CheY: 



T 



Y 



T P x + B UF^ T U x + B 



PF 



(13) 
(14) 



While CheB-P dephosphorylates spontaneously, 
the CheY-P hydrolysis is enhanced by the phos- 
phatase CheZ, an effect that is included in the high 
hydrolysis rate k HY for CheY-P (Lukat, 1991): 



B 



Y r - 

PF k 1 



Y, 



D 



UF 



The kinetic equations for these reactions are: 

rP] 4 



djY 1 
dt 



= Y, k n Y [TnW U \-k HY [Y P \, 



(15) 
(16) 

(17) 



n=0 



d\B 



PF] 



dt 



Y^k PB \T P ][B UF ]-k HB [B p % (18) 



n=0 



d[T P 
dt 



= k p [T^} - [T n p ][Y u ] - k PB [T p ][B^) 



J p - T p 



(19) 



where [Y u ] = [Y T ] - [Y p ], [B UF ] = [B F ] - [B PF ] 
and [TV] = [T n ] - [T p ]. J p is the net phospho- 
rylated receptor flux between methylation level n 



and (n + 1), given similarly as for J n in Eq. |Tl| 
with the free receptor concentration [T F ] replaced 
by the phosphorylated free receptor concentration 
[T„ ]. In all the above equations, the depen- 
dence on A is omitted, so the autophosphorylation 
rate and the phosphate transfer rates should all 
be considered as the rate for the combined recep- 
tor state (ligand occupied and unoccupied): k p = 
Lnk no + (1 — L n )k nv , k n — L n k no +{1 — L n jk nv , 
k PB = L n k PB + (1 — L n )k PB . It is also assumed 
that only CheB-P can bind with the receptors, 
which leads to the equation relating different sub- 
species of CheB: 



[B T ] = [B p ] + [B F ] - [B 



PF] 



(20) 



To describe the kinetics of the signal transduction 
pathway in full, we need to consider the interactions 
among the concentrations of all the 65 states for the 4 
chemical species: 60 receptor states = 2 ligand binding 
states x 5 methylation states x 3 enzyme binding states 
x 2 phosphorylation states, 1 free CheR state, 2 free 
CheB states and 2 CheY states. Using the fact that lig- 
and binding kinetics is fast and the enzymatic reactions 
are governed by Michaelis-Menten kinetics, the number 
of independent receptor concentrations is reduced from 
60 to just 10, consisting of the 5 free methylation states 
and the 5 phosphorylation states. Now, the whole system 
is described by kinetic equations Eq. Ofland Eq. [17 IK] 
plus conservation equations given by EqT0j9j and Eq. 2^7 

Concentration of the phosphorylated CheY ([Y^), 
which determines the tumbling frequency of bacteria, can 
be considered as the output of the whole chemotaxis sig- 
nal transduction pathway. In the next section, we study 
how the steady state concentration of CheY-P depends 
on the external ligand concentration [L], in particular, 
we derive a set of conditions for [Y p ] to be independent 
of [L], i.e., perfect adaptation. 



III. CONDITIONS FOR PERFECT 
ADAPTATION 

All the concentrations in our model fall naturally into 
two categories: the local variables defined for one partic- 
ular methylation level, such as [T n ] , the concentration of 
receptors with n methyl groups, and the global variables, 
such as [R F ], the concentration of the free CheR. The 
system adapts by adjusting the local variables with the 
ligand concentration, e. g., the steady-state values of [T n ] 
varies with [L]. However, perfect adaptation is achieved 
when the equilibrium value of [Y p ], a global variable, is 
independent of the ligand concentration (Othmer, 1998). 
This is generally not possible because the global variables 
are coupled with the local ones. One goal of this paper 
is to discover the conditions under which [Y p ] becomes 
independent of L. 

The strategy in obtaining the perfect adaptation con- 
ditions is to consider only global equations, such as the 
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conservation equations of the chemical species (e.g., Eq. 
fj], H and 20) and the steady-state equations of global vari- 
ables (e.g., Eq. [It] and [l8]), which do not depend on any 
one specific methylation level. In these global equations, 
there is no explicit dependence on ligand concentration, 

and composite variables, such as J2n=o^W^ m -^q. @, 
enter as weighted sums of the methylation level specific 
receptor concentrations. Another kind of global equation 
can be constructed by summing steady-state equations at 
all methylation levels (e.g., Eq. ||, [ll] and |l9|). The price 
to pay for such global equations is the introduction of 
new composite variables. However, if the reaction rates 
involved in different reactions are related in certain ways, 
the same composite variables appear in different global 
equations so that there are enough global equations to 
determine all the independent global and composite vari- 
ables. In other words, if certain conditions between reac- 
tion rates are satisfied, the steady-state concentrations of 
all the global and composite variables including [Y p ] can 
be independent of the ligand concentration, i.e., perfect 
adaptation. 

We leave the detailed derivation for the perfect adap- 
tation conditions to the appendix. In the following, we 
list these conditions, discuss their meaning and compare 
them with those found in previous works (Barkai, 1997; 
Yi, 2000). The perfect adaptation conditions can be 
grouped for each of the three pathway processes: condi- 
tion [l] is for the ligand binding and unbinding, conditions 
are required for the methylation process and condi- 
tions 1^-^ are related to the phosphorylation process: 

1. The time scale for ligand binding is much shorter 
than the methylation and phosphorylation time 
scale. This condition allows us to neglect ligand 
binding/unbinding kinetics. 

2. The association rates between the receptor and the 
methylation/demethylation enzymes, CheR and 
CheB-P, are linearly related to the activity of the 
receptor and are zero for n = 4 and n = 0, re- 
spectively: K^ x oc P±\ — P n \ and oc 
P n \ — Pq\. The dissociation rates of the enzyme 
receptor bound states are independent of A. 

3. The receptor activities of the non-methylated and 
the maximally methylated receptors are indepen- 
dent of A: P 0v = P 0o , P 4l , = P 4o . 

4. The ratios between the CheR catalytic rate (k^) 
and the CheB-P catalytic rate of the next methy- 
lation level (fc^ +1 ) are the same for all methylation 
states n: k^ +1 /k^ = const.. 

5. The phosphate transfer rates from Che A to 
CheB or CheY are proportional to CheA auto- 
phosphorylation rate: k p ^ oc P nX , k p ^ oc P nX . 



6. The explicit dependence on [T, F ] distribution can 
be removed from the expression 



[R F ] [B } 



K R 



K B 



n=0 



(21) 



This condition can only be strictly satisfied when 

[R F ] _ [B PF ] 

Condition |l| is necessary to decouple the ligand bind- 
ing process from the rest of the reactions. This is veri- 
fied experimentally and assumed in all the previous mod- 
els (Barkai, 1997; Morton-Firth, 1998; Spiro, 1997; Yi, 
2000). 

Condition 2 for the methylation process requires that 
the CheR and CheB methylation/demethylation rates 
depend linearly on the receptor's auto-phosphorylation 
rate (activity) . This is a generalization of the key ingre- 
dient for perfect adaptation found in Barkai and Leibler's 
work (Barkai, 1997). In the special case of Pn\ — 1 and 
Pox = 0, condition [2] means that CheB-P only bind to ac- 
tive receptors and CheR only bind to inactive receptors, 
the latter is missed in the original work of BL and later 
found to be necessary for perfect adaptation in (Morton- 
Firth, 1999) through a direct numerical simulation of the 
full system. 

The requirement in condition 3 that Pq\ and P±\ be 
independent of A is needed so that both the ligand-bound 
and vacant receptors have the same range of activity. 
This requirement for perfect adaptation is necessary in 
case the extreme methylation states n = or n = 4 
become populated with receptors. 

Condition 4 was first pointed out in (Yi, 2000), it is a 
more general form of the assumption that both fe„ and 
k^ are independent of n made in the original BL model. 
The justification of this condition may be related to a 
common evolutionary origin of CheR and CheB, result- 
ing in a similar anchoring position to the receptor for 
CheR methylating site n and CheB-P demethylating site 
n + 1 (Shapiro, 1994; Shapiro, 1995; Djordjevic, 1998; 
Barnakov, 1999). 

Condition 5 for the phosphorylation process is very 
similar to condition 2, in the sense that the phosphate 
transfer rates of the receptors have to be linearly related 
to their activity. This condition was not discovered be- 
fore because the phosphorylation process was neglected 
in previous works (Barkai, 1997; Yi, 2000). 

Condition 6 can only be satisfied exactly when one 
tunes the parameters such that the pre-factor in front of 
the sum in Eq. ^is zero. This condition was overlooked 
by most of the previous studies because the activities 
of the CheR or CheB-P bound receptors were neglected. 
However, in equilibrium, the population of enzyme bound 
receptors can be as high as 30% (Morton-Firth, 1999). 

By imposing all the conditions above, the steady state 
concentrations of the global variables will be independent 
of the ligand concentration, and are determined by 15 
parameters: the 4 total concentrations of table |, and 11 
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TABLE III: System parameters, numerical values from 
(Morton-Firth, 1999). 



Symbol Description Value 







1 2 


3 


4 




Relative activity of T n \ y 


0.125 0.5 


0.874 


1 







0.017 0.125 


0.5 


1 



K H CheR Michaelis constant 0.364 

K B CheB Michaelis constant 1.405 /j,M 

k R CheR catalytic constant 0.819 s _1 

k B CheB catalytic constant 0.155 s^ 1 

k p Che A autophosphorylation rate 15.5 s _1 

k PY CheA— >CheY phosphorus transfer rate 5/iM _1 s _1 

k PB CheA^CheB phosphorus transfer rate 5 /iM _1 s _1 

k HY CheY dephosphorylation rate 14.15 s _1 

k HB CheB dephosphorylation rate 0.35 s _1 



TABLE IV: Protein concentrations (in /iM) at different ligand 
occupancy rates L for the reference parameters. 



Species 


L = 


L = 0.5 


L = 1 


[To] 


0.028 


0.025 


0.002 


Pi] 


0.605 


0.316 


0.089 


[T 2 ] 


1.104 


0.923 


0.637 


[Ta] 


0.637 


0.947 


1.159 


[T A \ 


0.072 


0.289 


0.613 


[T A ] 


1.257 


1.250 


1.274 


[T P ] 


0.202 


0.201 


0.204 


[R F ] 


0.050 


0.050 


0.050 


[B F ] 


1.603 


1.602 


1.603 


[B PT ] 


1.858 


1.857 


1.860 


[B PF ] 


1.191 


1.190 


1.193 


\Y P ] 


1.200 


1.196 


1.209 



reaction rates of table [II, including P4 and Po, but not 
the relative activity values for the rest of the methylation 
states. However, for real biological system, these condi- 
tions for perfect adaptation may not be strictly satisfied. 
In order to understand bacteria's ability in adapting ac- 
curately under different internal and external conditions, 
i. e., robustness, we need to evaluate the effect of violat- 
ing these perfect adaptation conditions. 



IV. EFFECTS OF VIOLATING THE PERFECT 
ADAPTATION CONDITIONS 

Since it is not feasible to explore the whole parameter 
space, we choose to mostly perturb around the param- 
eter values that have been used in previous studies. To 
this end, we take most of our parameters from (Morton- 
Firth, 1999), which arc listed here in Table II and Table 
III. Hereafter we refer to this set of parameters as the 
reference parameters. Assuming ligand occupancy rate 
L n = L is independent of n, the steady state receptor 
distributions in different methylation states for different 
ligand occupancy rates L is shown in Fig. || A for the ref- 
erence parameters. In Fig. |]P, the population- weighted 
average receptor activities P n (L) = P no L + P nv (l — L) 
for methylation level n £ [0,4] is also shown. As is 
clear from Fig. Q when ligand (attractant) occupancy 
rate increases, the average receptor activity P n (L) de- 
creases for each methylation level n, and the system 
adapts by shifting the receptor population towards higher 
methylation states in achieving constant total activity 
[T A ] = Yl=o P n{ L )[ T n]- The steady state concentra- 
tions of all the other relevant concentrations at 3 dif- 
ferent ligand occupancy fractions are given in Table IV 
for the reference parameters, the small changes in [Y p ] 
at different ligand concentrations are caused by viola- 
tion of conditions 5 and 6 in the reference model used in 
(Morton-Firth, 1999) as we explain later in the section 



"Violating condition 5" . 

We have also constructed another model by modifying 
some of the reference parameters so that all the perfect 
conditions are satisfied. The results of perturbing this 
new model are essentially the same as for the reference 
model, mainly because the adaptation error in the refer- 
ence model is very small (< 1%). While this new model 
is mathematically more rigorous for isolating different er- 
ror sources, the reference model has the advantage that it 
is motivated biologically (from experiments or common 
sense), and therefore serves as a better starting point in 
exploring the parameter regions that are more likely to 
be biologically relevant. To make sure violation of condi- 
tions 5 and 6 in the reference model does not contaminate 
the effect of other conditions too much, we have always 
checked the error with and without violating the condi- 
tion in consideration, and made sure most of the error 
does come from violating the perfect condition we study. 

Since ligand binding is much faster than other relevant 
processes of the system, we do not consider the unreal- 
istic situation of violating condition 1. In the follow- 
ing, we study the effects of breaking the other 5 perfect 
adaptation conditions. Our goal is to understand the 
general reason behind the robustness of the system with 
respect to breaking each perfect adaptation condition. 
Even though we primarily perturb the system around 
the reference parameters, we also explore other parame- 
ter regions, especially when the reference model becomes 
insensitive to violation of a given condition. This strat- 
egy allows us to gain the general understanding of where 
in the parameter space a given perfect adaptation condi- 
tion becomes important and the reason behind it. 



A. Violation of condition 2 

Condition 2 requires that the methyla- 
tion/demethylation enzyme binding rates to a receptor 
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1 2 3 

Methyl groups 



1 2 3 

Methyl Groups 



FIG. 2: (A) Distribution of receptors in different methylation states at different ligand occupancy fractions L for the reference 
parameters. The total activity of the system being [T A ] — 0.5[T T ]. (B) The population-weighted average receptor activity 
Pn(L) for different methylation level n € [0,4] at different fractional ligand occupancy rates L. 



depend linearly on the activity of the receptor. For 
the reference parameters, where Pq — and P4 = 1, 
condition 2 simply means that CheR only binds to 
inactive receptors and CheB-P only binds to active 
receptors. The simplest way in violating condition 2 is 
to allow CheR bind to active receptor or CheB-P bind 
to inactive ones , which can be formally expressed as: 



K* 1 = K R V(l - P nX + Or), 



K 



= K' 



(22) 



where a r > and > are the measures of violating 
condition 2; b r and bb are normalization factors tuned 
with respect to a r and at, to keep the total activity of 
the system constant at a given ligand occupancy rate 
(L=0.5) for comparison purpose. a r = and a& = cor- 
responds to condition 2 being satisfied; a r — > 00 (with 
a r b r = const.) or a& — > 00 (with abbb = const.) respec- 
tively corresponds to CheR or CheB-P binding to all re- 
ceptors equally. 

In Fig. ||, we show the steady-state concentration of 
CheY-P versus the ligand occupancy rate L for vari- 
ous values of a r and at,. Even for the extreme cases 
of a r — 00 or ab = 00, respectively corresponding to 
CheR or CheB-P binding to both active and inactive re- 
ceptors equally, the deviation from perfect adaptation is 
only ~ 10% — 15%. Intuitively, the reason for the near 
perfect adaptation is that the control of the system's 
total activity can be carried out by either the methy- 
lation (CheR) or demethylation (CheB-P) process, pro- 
vided that at least one of the enzymes' binding rates 
is strongly correlated with the receptor activity. If the 
receptor binding rates of both enzymes become indepen- 
dent of the receptor's activity, i. e., both a r = 00 and 
a,b = 00, the system is only controlled through the weak 



effect of CheB phosphorylation and does not adapt very 
well. 

Specifcally, condition 2 requires that CheR does not 
bind to the fully methylated receptors (n=4), and ChcB- 
P does not bind to the unmethylated receptors (n=0). 
Therefore, the quantitative effects of breaking condition 
2 (as in Eq. p2| ) depends on the receptor concentration at 
the fully methylated state [T4] or the unmethylated states 
[To] (see appendix for details). Both [To] and [T 4 ] are 
relatively small for the reference parameters with [T4] > 
[To] (see Fig. [2]), which explains the qualitative features 
in Fig. |. The effect of — > 00 only becomes noticeable 
because To] is not too small for — > 00. 



B. Violation of condition 3 

Since adaptation for bacterial chemotaxis relies on bal- 
ancing the effect of ligand binding on the receptor's ac- 
tivity with that of the methylation of the receptor, a 
necessary condition for perfect adaptation is for both lig- 
and bound and vacant receptors to have the same range 
of activity, i.e., condition 3. For the reference parame- 
ters, condition 3 is obeyed by having: Pq v = Pq = 0; 
Pi V = P40 = 1 • Without changing the monotonic depen- 
dence of the receptor activity on their methylation level, 
we can break condition 3 at n = by increasing Pq v from 
to |; or at n — 4 by decreasing P± a from 1 to |. The 
enzyme binding rates are adjusted accordingly in keep- 
ing condition 2 satisfied. The effects are shown in Fig. |j. 
The system is insensitive to the opening of the activity 
gap AP = Pov — Poo at n = 0, because the receptor 
population is small at n = even at L = 0. For the 
same opening of activity gap AP4 = P iv — P 4o at n = 4, 
the adaptation error is 6%. In particular, the system has 
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FIG. 3: The steady state [V ] concentration versus ligand binding rate L for different ways of breaking condition 2: (A) CheB-P 
binds with active receptors only (a& = 0), CheR is allowed to bind with active receptor with varying strength a r — 0, .2, 1, 5, oo, 
where a r = oo corresponds to CheR binds to all receptor indiscriminately. (B) CheR binds with inactive receptors only 
(a r = 0), CheB-P is allowed to bind with inactive receptor with varying strength ab = 0, .2, 1, 5, oo, where ab — oo corresponds 
to CheB-P binds to all receptor indiscriminately. 



lower CheY-P concentration at higher ligand occupancy 
rate L, because the receptor population shifts towards 
higher methylation levels at larger L, and the effect of 
methylation is not large enough to cancel the decrease 
of activity caused by ligand binding. Quantitatively, the 
adaptation error increases with the activity gap; e.g., it 
reaches 25% when we lower P^ Q further to 0.5. 



C. Violation of condition 4 

The methylation and demethylation catalytic rates k3 
and fc„ can depend on methylation level n. From Eq. |llj, 
the steady-state properties of the system only depend on 
the ratios: r n = k^/k^_ 1 for n — 1,2,3,4. Condition 
4 for perfect adaptation requires that r n be a constant 
independent of n, a kind of "detailed balance" condition. 
Indeed, if we change k^ and k^ while keeping r n con- 
stant, the system adapts perfectly. However, when we 
make r n depend on n, perfect adaptation is lost. In Fig. 
||, we show the effects of increasing one r n by a factor of 
2 while keeping the rest r n unchanged at their reference 
value for n = 1,2,3,4 respectively. The quantitative de- 
viation from perfect adaptation depends on n, with the 
largest deviation of ~ 25% occurring at n = 2 , possibly 
because the receptors are highly populated at n = 2 for 
the reference parameters. 



D. Violating condition ^ 

Condition 5 requires that the phosphate transfer rates 
of a receptor is proportional to its auto-phosphorylation 




FIG. 4: The steady state [Y p ] concentration versus ligand 
occupancy fraction L for different ways of breaking condition 
3 by opening the activity gap at n — 0: Po v = Piv = g (long 
dashed line) or at n — 4: P^ = Ps v = ~ (short dashed line), 
the solid line is for the reference parameters. The two inserts 
illustrate the opening of the activity gap at n = and n = 4 
respectively. 



rate, a kind of compatibility condition. The simplest 
way to break condition || is to set the phosphate transfer 
rates to be a constant independent of both the ligand 
binding and the methylation level of the receptor. This 
assumption is also made in (Morton-Firth, 1998; Morton- 
Firth, 1999). 
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FIG. 5: Steady-state CheY-P concentration [Y p ] versus lig- 
and occupancy rate L for breaking condition 4. The results 
are obtained by increasing one of the 4 ratios of catalytic rates 
(see text for definition) by a factor of 2: r n = 2r re / = 0.38, 
while keeping the other 3 ratios unchanged at the reference 
value of 0.19. The 4 curves correspond to n = 1,2,3,4 re- 
spectively. 



For the reference parameters, the steady-state [Y p ] 
change by less than 1% over the whole range of ligand 
occupancy as shown in Fig. |^ (curve a), indicating the 
insensitivity of the system's perfect adaptation with re- 
spect to this particular choice of breaking condition In 
the following, we explain the system's near perfect adap- 
tation by the existence of approximate global equations. 

In deriving condition ^, a global equation is formed by 
summing Eq. |l9| o ver all methylation levels and replacing 
[T%] by [T n ] -[T* ], which leads to the formation of 4 com- 
posite variables: G = En=0 fe ™[ T «]> G l = J2i=0 k n [ T n], 
G2 = £*=o*£ r Kf] and G 3 = ELo^K}- Condi- 
tion [5| is needed to make G±, G2 and G3 proportional 
to each other, so that the total number of global equa- 
tions is enough to solve for all the independent global and 
composite variables (see section "Conditions for perfect 
adaptations" and Appendix for details). When condition 
U is broken by setting k PB and k PY to be constant, G2 
and G3 are still proportional to each other, but they are 
now different from G\ , the total number of global equa- 
tions are now not enough in solving for all the global vari- 
ables, and local equations have to be used. This leads to 
all the global variable depend on ligand concentration, 
i.e., non-perfect adaptation. However, because the con- 
centration of (un-phosphorylatcd) CheY is much larger 
than the receptor concentrations, the phosphorylated re- 
ceptor concentration [Tjf] is small compared with the 
total receptor concentration [T n ], due to efficient phos- 
phate transfer from CheA to CheY and the subsequent 
high CheY-P dephosphorylation rate. As a result, G\ 
is negligible relative to Go, leading to an approximate 



FIG. 6: Relative steady-state CheY-P concentrations 
rF p ]/[Y p ]i = o versus ligand occupancy rate L when condition 
is violated; the adaptation error depends on the parameters 
of the system, 4 cases are studied here using parameters with 
increasing degrees of deviation from their reference values (see 
text for detail). The parameters used are: (1) For curve a, 
reference parameter values; (2) For curve b, same as for curve 
a, except [Y T ] = [Y T ]R C f/20; (3) For curve c, same as for 
curve 6, except k R — fcj^ of /50; (4) For curve d, same as for 
curve c, except P„i = 0.25 and P V 2 = 0.6. 



global equation with the same degree of reduction in in- 
dependent composite variables and eventually the near 
perfect adaption observed in Fig. || A. 

However, reducing CheY concentration alone does not 
change too much the system's ability in perfect adapta- 
tion, as shown in Fig. (curve b). At low CheY concen- 
tration, the phosphate group of CheA-P goes to CheB. 
Because of the slow dephosphorylation rate of CheB-P, 
most of the CheB become phosphorylated in steady state, 
essentially decoupling the phosphorylation process from 
the adaptation process. The adaptation of the system 
therefore becomes insensitive to the phosphorylation re- 
lated condition 5. 

To amplify the effect of violating condition 5, we re- 
duce the overall activity to [T A ] = 0.014[T T ] at L = 

r making k R = 0.02/c|? ef . The result is shown in Fig. 
(curve c). The adaptation accuracy can also depend 
on other parameters, such as the receptor activity P n \. 
In Fig. |6j (curve d), we show that a slight change in 
receptor activity leads to higher deviation from perfect 
adaptation. 



E. Violating condition 

The total receptor activity [T A ](= £* =0 -P„[T„]) is di- 
rectly related to the final production of CheY-P. How- 
ever, only part of [T A ] can be expressed in terms of other 
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composite variables related to receptor population, i. e., 
the total free receptor concentration [T F ] = Xm=nP?f ] 
and the total activity due to free receptors [T" 4 ^] = 

X^n=o-fn[^rf]- ^ nas an ex tra term £ coming from the 
activity of the enzyme (CheR or CheB-P) bound recep- 
tors (see appendix for details), which is proportional to 

= Eto^M with a pre-factor (-Igl + 1^1) 
(see Eq. |2l]) . Condition || is required to eliminate this 
extra global variable £' by setting the pre-factor to zero. 

The effect of breaking condition || can be small, be- 
cause as [R F ] deviates from its perfect adaption value 
[R F ]Ada P , so does [B PF ] with the same trend, leading 
to small changes of the pre-factor in £. Also, part of £' 
can be approximated by a linear combination of [T F ] and 
[T ] , depending on the activity levels of different recep- 
tors P n \- Finally, for higher total activity, the relative 
effect of £ will be small. For the reference parameters, 
the accuracy of adaptation is better than 98% for 4-fold 
change of CheR concentration from its perfect adapta- 
tion value, as shown in Fig. ^ A. The adaptation accu- 
racy decreases as we lower the total activity by decreas- 
ing methylation rate as shown in Fig. [?] B. Finally, 
when we increase the activity differences between the lig- 
and bound and the vacant receptors by setting: P no = 
(n = 0, 1, 2, 3), P 4o = 1; P nv = 1 (n = 1, 2, 3, 4), P 0v = 0, 
the same change in [R T ] can cause more than a 50% error 
in adaptation, as shown in Fig. ffl C. 



V. COMPARISON TO STOCHASTIC 
SIMULATION AND EXPERIMENTS 

The results from the previous sections can be compared 
with both the discrete stochastic numerical simulation 
and real experiments. We use the reference parameters 
for all the comparison studies. 



A. Comparison to stochastic simulation 

Stochsim (Morton-Firth, 1998) is a general purpose 
stochastic simulator for chemical reactions. For our 
study, the volume of Stochsim simulation is set to 
be 1.4" 15 L, and the number of molecules is therefore 
843 x concentration (in /xM). 

In Fig. H A, we show the Stochsim simulation result 
for the reference parameters, which agrees well with the 
results from simulating our continuum equations with the 
same parameters. In Fig. || B, we show the Stochsim 
simulation result for the parameters used in Fig. ^ C 
with [R T ] — 4:[R T ]Adap, where perfect adaptation is lost 
because of violation of condition 6. As predicted from 
our deterministic model, after sudden changes of ligand 
occupancy rate L, [Y p ] does not always return to its 
pre-stimulus level; in fact, the maximum error (~ 50%) 
is observed when L = 0.2, consistent with Fig. [7] C. 

For most of the results shown in this paper, we have 



compared with the results from stochastic simulation us- 
ing Stochsim (data not shown). Overall, the averaged 
behaviors of Stochsim simulations are consistent with our 
continuum model, which is interesting given the non- 
linear nature of the chemical kinetics. Further work is 
needed in characterizing the fluctuation of the individual 
Stochsim simulations, and compare them with the fluc- 
tuations in behavior among different individual bacteria 
(Morton-Firth, 1998). 



B. Comparison with experiment 

In a recent experimental study by Alon et al (Alon, 
1999), mutant bacteria lacking a certain chemotaxis pro- 
tein, such as CheR, CheB, CheY or CheZ, are used, and 
the missing protein is reintroduced in a controlled fashion 
through a plasmid inserted into the mutant bacteria cells. 
This technique allowed these authors to study the effect 
of various enzyme concentration changes on the chemo- 
taxis behavior of the bacteria. Specifically, the tumbling 
frequency of the bacteria is measured through a sudden 
increase of ligand concentration, which effectively corre- 
sponds to a sudden change of ligand occupancy rate from 
L = to L = 1. 

In Fig. H A, we show the adaptation precision as 
the ratio between phosphorylated CheY level before and 
after the stimulus for various CheR and CheB concen- 
trations. For CheR concentration change of up to 50 
fold with respect to the reference value, the adaptation 
error is < 3%, somewhat smaller than the experimen- 
tally measured adaptation error cited in (Alon, 1999). If 
[B T ], instead of [R T ], is changed, the adaptation error 
would be much bigger, as shown in Fig. ^| A. This is 
the case because for large values of [B T ], the low activ- 
ity and the large values of [i? pi? ] make the violation of 
condition 5 and 6 more significant. This could explain 
the larger (1.09) adaptation precision reported in (Alon, 
1999) when [B T ] expression is 12 times that of the wild 
type values. Since we define adaptation accuracy based 
on CheY-P concentration, the quantitative difference be- 
tween the adaptation error observed in (Alon, 1999) and 
those of our model could be explained by the signal am- 
plification at the motor level (Cluzel, 2000). 

The relaxation time of the system after a sudden 
change in ligand concentration can be determined by 
direct simulation of the full kinetic equation or by lin- 
earizing the methylation/demethylation kinetic equa- 
tions around the steady state. The dependence of both 
the steady-state tumbling frequency and the linear re- 
laxation time on CheR concentration [R T ] is shown in 
Fig. | B. They agree qualitatively with the steady state 
tumbling frequency and the relaxation time measured in 
(Alon, 1999), as depicted in Fig. 2 B of their paper, 
although direct quantitative comparison is not possible 
due to different definitions of relaxation time and lack 
of detailed understanding on how CheY-P regulates the 
motor. 
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FIG. 7: Relative steady-state CheY-P concentrations [Y ]/[Y ]l=o versus ligand occupancy rate for different CheR concen- 



trations (with condition 5 satisfied), which are varied with respect to the perfect adaptation value [R ]Ada P - (A) Reference 
parameters are used except the different values of [R T ] listed in the figure, and [R T ]Ada P = 2.63[i? T ] re /; (B) k R = 0.1fc^,j is 
chosen in reducing the total activity, where adaptation is less accurate, and [R T ]Ada P = 5.26[i? T ] re /; (C) Same parameters as 
in B except that the activity difference between ligand bound and vacant receptors are set to be maximum (see text), and 
[R T ]Ada P = 5.35[7? T ] re /. 



VI. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have studied a theoretical model de- 
scribing the full chemotaxis signal transduction pathway. 
Through systematic analysis of the steady state proper- 
ties of the model, we derive a complete set of conditions 
for the system to adapt exactly. Some of the conditions 
are generalizations of the ones discovered before, but oth- 
ers, in particular, the conditions related to the phospho- 
rylation part of the pathway, are discovered for the first 
time here. It is quite remarkable that perfect adaptation 
can be achieved for arbitrary ligand concentration with a 
small set of conditions, far less than the number of vari- 
ables and the number of reaction rate constants in the 



problem. 

The (intrinsic) state of a receptor can be described 
by its ligand binding status (A) and mcthylation level 
(n). The (external) properties of the receptor com- 
plex include its abilities to interact with the methyla- 
tion/demethylation enzymes, to undergo autophosphory- 
lation, and to transfer its own phosphate group to CheY 
or CheB, all of which depends on the (internal) state of 
the receptor characterized by n and A. Perfect adapta- 
tion requires these three properties of the receptor com- 
plex to be correlated with each other in a linear fashion 
for any given receptor state {n\} (condition 2 and 5). 
Available experimental data that addresses the validity 
of such connections has been discussed extensively in (Yi, 
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FIG. 8: Dynamics of [Y p ] from Stochsim simulation with ligand occupancy rates L changing from 0^1^ 0.2 at 50 and 
250 seconds when the parameters are set to: (A) the reference values, (B) same as in Fig. Q C with [R T ] = 4{R T ]Ada P - The 
solid lines are results from simulations of our deterministic equations, the dotted lines are fits to the Stochsim data with an 
exponential decaying function to obtain the relaxation time. 
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FIG. 9: The response to a sudden increase of ligand concentration determined from the continuum model. (A) The steady- 
state CheY-P concentration ratios before and after the stimulus [Y p ]l=o/[Y f ]l=i for different fold changes of CheR or CheB 
concentrations; (B) Steady state CheY-P concentration and the linear relaxation time upon sudden change of ligand occupancy 
rate (from to 1) versus different CheR concentrations. 



2000). Even though the evidence for such connections is 
not well established and the correlation may not be lin- 
ear, it is conceivable that a high degree of correlation 
exists among these three properties of the receptor, be- 
cause they are determined by the same conformational 
change of the receptor protein complex for a given recep- 
tor state {n\}. 

Since most of the perfect adaptation conditions are 
relations between different reaction rates, the system's 
ability to adapt accurately can be considered "robust" in 



the sense that the perfect adaptation is independent of 
concentrations of any specific chemotaxis protein, which 
can fluctuate between different individual cells and at dif- 
ferent stage of the cell development. Only one of the per- 
fect adaptation conditions requires the fine tuning of the 
methylation enzyme concentrations (condition 6). Be- 
cause of this condition, in the strict mathematical sense, 
the perfect adaptation of the system can only be achieved 
via fine tuning of a parameter, and therefore cannot be 
considered robust. However, as we have shown in this 
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paper, the effect of violating this condition can be rather 
small, especially at the reference parameters. 

The discovery of the perfect adaptation conditions pro- 
vides an invaluable starting point in exploring the param- 
eter space. We evaluate the sensitivity of the system's 
perfect adaptation ability by perturbing the perfect adap- 
tation conditions. We find that the system can adapt 
near perfectly even in the absence of some of the per- 
fect adaptation conditions. In finding the perfect adap- 
tation conditions, we focus on studying equations which 
do not depend on any individual mcthylation levels, these 
global equations are obtained by either conservation laws 
or summing steady-state equations over different methy- 
lation levels. The same approach is also useful in under- 
standing the near perfect adaptation when the perfect 
adaptation conditions are violated. Technically, we can 
explain the near perfect adaptation by the existence of 
approximate global equations replacing the ones lost due 
to the violation of perfect adaptation conditions. Bio- 
logically, these approximate global equations are caused 
by various intrinsic properties of the system, such as sep- 
aration of scales in protein concentrations and reaction 
rates, or specific properties of the receptor distribution 
in different methylation states. Since real biological sys- 
tems are not likely to satisfy all the perfect adaptation 
conditions exactly, the abundance of such near perfect 
adaptation regions in the parameter space strongly lim- 
its the range of activity variation and is probably respon- 
sible for the robustness of the system's ability to adapt 
almost perfectly. 

Through systematic study of the system's behavior 
when different perfect adaptation conditions are violated, 
we have also identified parameter regions where signifi- 
cant deviation from perfect adaptation occurs. This may 
provide possible explanations to bacterial chemotaxis re- 
sponses that does not adapt accurately, such as the serine 
response as reported in (Berg, 1972), and constitute con- 
crete predictions that can be experimentally verified. 

Aside from perfect adaptation, another challenge for 
modelling bacterial chemotaxis is to understand the large 
signal amplification from ligand concentration change to 
the change in bacterium flagella rotation bias. To di- 
rectly compare between experiments and simulation, de- 
tailed information between CheY-P concentration and 
the motor rotation bias is needed. Recently, the con- 
nection between CheY-P level and the motor activity 
was investigated in (Scharf, 1998; Alon, 1998; Cluzel, 
2000). In (Cluzel, 2000), where rotation bias of single 
bacterium was measured for different [Y p ] concentra- 
tions, it was shown that the motor bias for individual bac- 
terium should be fitted by a Hill function with a large Hill 
coefficient (~ 10). This highly nonlinear function may 
explain the advantage of perfect adaptation in amplify- 
ing the gain, and also the nonlinear dependence of Bqcw, 
the CCW rotation bias, on changes in ligand occupancy 
as found in (Jasuja, 1999). However, quantitatively, from 
(Cluzel, 2000), the maximum signal amplification from 
change in [Y p ] to the tumbling frequency is measured to 



be: dB ccw /dln[Y p ] w 2.2. With the reference parame- 
ters in our model, this leads to a total signal amplification 
of dB C cw/dln[Y p ] x d\n[Y p ]/dL w 2.2 x 0.65 « 1.43, 
which is still much too small as compared with the total 
signal amplification measured in experiments, e.g., ~ 30 
as reported in (Jasuja, 1999). 

The gain of the system could come from receptor clus- 
tering as suggested in (Bray, 1998). However, to reconcile 
the existence of high gain and the wide dynamic range of 
response, it is highly desirable to have high gain for the 
signal transduction pathway itself. One of the interest- 
ing findings of our study is that if the system satisfies all 
the perfect adaptation conditions, the steady state ac- 
tivity of the system is independent of the exact values 
of the receptor activity P n \ for n e [1,3]. On the other 
hand, the response of the system, defined here as the 
difference of CheY-P concentrations between its extreme 
value after the stimulus and its original value before the 
stimulus, directly depends on the difference of receptor 
activity between ligand bound and ligand free receptors: 
AP n = P nv — P no . The higher these differences are, the 
higher the response will be. In order to have high re- 
sponse, it is favorable to increase AP n and to have lower 
total activity. Indeed, if we simply increase the activity 
difference between the ligand bound and vacant receptor, 
such as those used in Fig. c, the total amplification can 
be increased to: 2.2 x 1.7 = 3.74. Other changes, such as 
reducing the system's total activity, can enhance the gain 
much more, as noted also in (Barkai, 2001). A detailed 
study of the response of the system is outside the scope 
of this paper and will be reported in another communi- 
cation. 

Overall, the current model is capable of explaining the 
qualitative behaviors of the chemotaxis pathway related 
to adaptation, in particular, the robustness of the sys- 
tem's ability to adapt nearly perfectly. Much work is 
still needed to modify and enrich the model to under- 
stand the high sensitivity and wide dynamic range of the 
system (Sourjik, 2002). However, because adaptation 
and response occur with very different time scale and via 
largely different molecular processes, modification of the 
model in explaining the high response gain should not 
change the perfect adaptation conditions significantly. 
Indeed, it is not hard to show that even with receptor 
coupling added to the current model, the conditions we 
identified in this paper are still needed for the system to 
achieve perfect adaptation, the only change is that activ- 
ity of each receptor now depends also on its neighbors' 
activities (B. Mello and Y. Tu, manuscript in prepara- 
tion). We believe that, as long as the basic structure of 
the protein interaction network stays intact, the perfect 
adaptation conditions identified here will be mostly valid. 
These conditions not only offer explanation for adapta- 
tion accuracy and its robustness, furthermore, they serve 
as constrains for constructing quantitative models in un- 
derstanding other aspects of the bacterial chemotaxis. 
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Appendix 

In this section, we describe the detailed derivation of 
the perfect adaptation conditions listed in the section 
"Conditions for perfect adaptation" .As described there, 
the approach is to construct global equations using global 
and composite variables that do not depend on the recep- 
tor population in any one individual methylation state. 

First, we concentrate on the methylation related equa- 
tions. Eq. 0, Eq. |^ and summation of Eq. ^ over 
n € [0,4] gives 3 global equations. For the steady state, 
the methylation flux between different methylation states 
should be zero: 



J' n 



K R 



[b pf ][t f +1 ] 



0, (0<n<3). 

(23) 

Condition H can be used in factoring out the common n 
dependent factor from k p +1 and k B in J„, after which 
Eq. ^3] are summed over n 6 [0,3] to obtain a global 
equation. 

Using condition ||, the Michaelis constants can be 
expressed as K R X = K R /(P ix - P nA ) and K B X = 
K B /(P n \ - P x), where K R and K B are constants. If 
we further enforce condition 3, i.e., P± a = P^ v = P4 and 
Pq = Pq v = Pq, we can convert all the weighted sums 
of the individual receptor concentrations into two com- 
posite receptor concentrations [T F ] and [T AF ]. [T F ] = 
J2 n= o[T F ] is the total concentration of the free receptors; 
[T AF ] = En=o P »[ T n] is the total concentration of the 
active free receptors, where P n = (1 — L n )P nv + L n P no is 
the population-weighted average activity for a receptor 
with n methyl groups. Therefore, after applying condi- 
tions 2, 3 and 4, the 4 methylation related global equa- 
tions can be written as: 



[R T ] = [R F }(1 + P 4 
[B P } = [B PF }(1-P ( 
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-P a [T F ] + [T AF ]) 



K B 



0. 



(24) 
(25) 

(26) 
(27) 



If is a constant (i.e. 



.ouanu ^i.c, CheB-P binds equally to all 
receptors), condition 2 is violated. However, it is not 
hard to see that if the receptor population in the n = 
methylation state, [To], is small, we can still sum up the 
methylation balance equations to form a global equation. 
The same is true if K R X is a constant and [T 4 ] « 0. 

Next, we focus on the phosphorylation related equa- 
tions. Besides its importance in producing the final 



output of the signal transduction pathway CheY-P, the 
phosphorylation is also coupled back to the methyla- 
tion process through concentration \B PF ]. By writing 
ri A and using condition ||: k PY = k PY P n \ 
k PB P n \, the phosphorylation related global 
equations can be written as: 



Kx = k 
and fc£f 



[Y p ] 
[B PF ] 



„PY \rpPA 



[T 1 



k HY + k PY [T PA ' 
k PB [T PA ] 



T PA ] = 



k HB + k PB [T PA ] 
\T A ] 



[Y T ], 
[B F ], 



(28) 
(29) 



1 + ^([Y T ] - [Y p ]) + ^([B F ] - [B PF ]Y 



(30) 



over n £ [0, 4]. 
Al - the 



Eq. |3^ is obtained by summing Eq. |l9| 
There are two composite variables, [T^A] and [T"j m 

above equations. [T A ] = 5Z n=0 -Pnpra] is the total con- 
centration of active receptors, [T PA ] = ^n=o.AP n \\T p x ] 

X—v,o 

is phosphorylated active receptor concentrations. 

If the CheA phosphate transfer rates are independent 
of its ligand/methylation status, i.e., k PY = k PY and 
k PB = k PB , condition 5 is broken. A new composite 
variable [T p ] = ^ n=0 [T p } appears in the above equa- 
tions, replacing [T Pj4 ] in Eq. [23, Eq. ^ and part of 
Eq. g However, if [T PA ] <c[T j4 ], e.g., due to effi- 
cient phosphate transfer from CheA to CheY, [T Pj4 ] can 
be neglected, and there again will be only two composite 
variables in the phosphorylation related global equations, 
and therefore the system may still adapt near perfectly 
in absence of condition 5, as discussed in the section "Vi- 
olating condition 5" . 

The methylation and the phosphorylation global equa- 
tions communicate through various CheB concentrations. 
An extra equation is necessary to connect the concentra- 
tions of these different forms of the same proteins: 



[B T ] = [B p ] + [B F ] - [B 



PFi 



(31) 



Finally, by using Eq. ^|, we can write down the expres- 
sion for the total receptor activity of the system [T A ] that 
appears in Eq. p0|: 



= 1 



[R F ] [B PF ] 

[R F ] [B PF ] S 1 
" K R + K B 



\T 



J2Pn 2 [T F }- (32) 



n=0 



The above equation contains a new composite variable 
£' = X)n=o^n 2 [^f]- Condition H is thus required to elim- 
inate this extra term. Part of £' can be expressed in terms 
of the other composite variables, such as [T F ] and [T AF ]. 
Therefore, the effect of violating condition 6 can not be 
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simply measured by the value of as we discussed in 
the section "Violating condition 6" . 

If all the conditions listed in Conditions for perfect 
adaptation are satisfied, we have nine global equations: 
Eqs. [j4[ - |32| , these 9 global equations contains 5 global 
variables: [R F ], [B p ], [B PF ], [B F ], [Y p ], and 4 com- 
posite variables: [T F ], [T AF ], [T A ], [T PA ]. Therefore, 
the steady state values of all the nine global or compos- 
ite variables, including [Y p ], will be independent of the 
ligand concentration and the system can achieve perfect 
adaptation. 

We are thankful to Jeremy Rice, Geofferey Grinstein 
and Gustavo Stolovitzky for helpful discussions and care- 
ful reading of the manuscript. B. Mello has a scholarship 
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